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ABSTRACT 



Aims. We report on archival Hubble Space Telescope (HST) observations of the globular cluster M71 (NGC 6838). 
Methods. These observations, covering the core of the globular cluster, were performed by the Advanced Camera for Surveys (ACS) 
and the Wide Field Planetary Camera 2 (WFPC2). Inside the half-mass radius (n, = l'.65) of M7I, we find 33 candidate optical 
counterparts to 25 out of 29 Chandra X-ray sources while outside the half-mass radius, 6 possible optical counterparts to 4 X-ray 
sources are found. 

Results. Based on the X-ray and optical properties of the identifications, we find 1 certain and 7 candidate cataclysmic variables 
(CVs). We also classify 2 and 12 X-ray sources as certain and potential chromospherically active binaries (ABs), respectively. The 
only star in the error circle of the known millisecond pulsar (MSP) is inconsistent with being the optical counterpart. 
Conclusions. The number of X-ray faint sources with Lx > 4x 10^" ergs s"' (0.5-6.0 keV) found in M71 is higher than extrapolations 
from other clusters on the basis of either collision frequency or mass. Since the core density of M7 1 is relatively low, we suggest that 
those CVs and ABs are primordial in origin. 

Key words, globular clusters: individual (M71, NGC 6838) 
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There are 158 Galactic globular clusters (GGCs) found in the 
halo of our galaxy and typically contain 10"^ - 10^ stars. They 
. are very old and dense star systems and tightly bound by grav- 

■ ity, which gives them their spherical shapes and relative high 
stellar density toward the center The dense stellar environment 
in globular clusters triggers various dynamical interactions, ie. 
exchanges in encounters with binaries, direct collisions, destruc- 

l tion of binaries, and tidal capture. These dynamical interactions 

■ not only can change the evolution of individual stars, but also can 
' produce tight binary sy stems (see, e.g., Ashman & Zepf 1998; 

IVerbunt & Lewin|[2004l for review). 

One of the most powerful ways to probe the binary content 
of globular clusters is by studying the X-ray source population. 
In the early 1970s, X-ray sources with luminosity greater than 
10^^ ergs s ' were first detected by using the Uhuru and OSO-7 
Observatories. Following the Einstein and ROSAT era, the num- 
ber of faint X-ray sources (Lx < 10^"^'^ ergs s"') was dramati- 
cally increased. Those bright X-ray sources have been identified 
with low-mass X-ray binaries (LMXBs; iGrindlav et al.l Il984l) 
while the identification of the weaker sources remained limited 
due to low photon statistics and insufficient spatial resolution. 
The launch of the Chandra X-ray Observatory ushered in a new 
age of studying the crowded centers of Galactic globular clusters 
with a far greater sensitivity and resolving power than ever be- 
fore (e.g., IGrindlav et ai] 12001 alfh) . With the aid of the Hubble 
Space Telescope (HST), many of these faint X-ray sources 



were identified as quiescent low-mass X-ray binaries (qLMXBs; 
in which a neutron star accretes matter from its companion 
at a low rate), cataclysmic variables (CVs; in which a white 
dwarf accretes from its low-mass companion), and millisecond 
pulsars (MSPs), as well as chromospherically active binaries 



(ABs; e.g., RS CVn and BY Dra systems) (e.g., IGrindlav et al 
200fbl:lPoolev et al.>2002; Edmonds et al.l2003airbl;lHeinke et al 



20051; iBassa et alj 12004 iKong et alj l200a iLugger et all 12007 
Bassa et al.ll2008l) . 



The globular cluster M71 (NGC 6838) lies close to the 
Galactic plane with Galactic longitude / = 56°74 and latitude 
b - -4°.56. Similarly to 47Tuc, it is a fairly metal-rich globular 
cluster with metallicity of [Fe/H] = -0.73. Due to its relatively 
small distance (d ~ 4kpc) to Earth and low central luminosity 
density (pc - lO"* "^ Lg/pc^), M71 is a good target for both opti- 
cal and X-ray observations. The core, half-mass, and tidal radii 
are = 0'63, = 1.'65, and =8f96, respectively. M71 shows 
no evidence for core collapse. Its moderate optical reddening 
Eb-v = 0.25 may be converted in to a nominal X -ray absorption 
column of A^H = 1.39x 10^'cm'^ ( iPredehl & Schmitt 1995) . The 
aforemention ed parameters related to M7 1 were obtained from 
iHrnisi ([19961 updated 200j|)- 

In this work we report on archival Chandra and HST ob- 
servations of the globular cluster M71. We have obtained a 
52.4-kilosecond Chandra observation of M71 taken with the 
Advanced CCD Imaging Spectrometer (ACIS), reaching the 
limiting X-ray luminosities of 1.5 x lO-'" ergs s"' and 7.9 x 
10^^ ergs s"' in the energy ranges of 0.3-8.0 and 0.5-2.5 keV, re- 
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spectively. In lElsner et alj (l2008l) . we reported the identification 
of 29 X-ray sources within the cluster half-mass radius, includ- 
ing the known milHsecond pulsar PSR J1953-I-1846A (M71A), 
and their X-ray properties, and found that 18 + 6 of these 29 
sources are likely to be associated with M71 from a radial distri- 
bution analysis. The present paper extends our study of the X-ray 
sources in M71 by using archival HST data to identify optical 
counterparts to the majority of M71's X-ray sources, improving 
our understanding of their nature. 

In §2 we discuss the Chandra X-ray observations and spec- 
tral analysis. HST observations, data reduction and analysis are 
described in §3. In §4 we present the source identification. A 
discussion and comparison with other globular clusters is given 
in §5. 

2. X-ray Observations 

lElsner et"aL I (12008 ) described the Chandra X-ray observations 
of M7 1 . We here note some relevant information for our anal- 
ysis and extend the spectral fitting in this paper to test alterna- 
tive m odels besides the power-laws considered by lElsner et"an 
(l2008h . Only seven of the 29 detected X-ray sources have suf- 
ficient counts (six sources with at least 50 source counts and 
one known MSP with 37.5 source countfl) to warrant a de- 
tailed spectral analysis. We used the CIAO tool dmextract to 
extract spectra of the brighter sources and the source-free back- 
ground regions near to those sources. Response files were con- 
structed by using the CIAO tool mkacisrmf and mkarf. The ex- 
tracted spectra were binned with at least 5 source counts per bin. 
Background-subtracted spectral modeling was performed with 
XSPEC using data in the energy band 0.3-8.0keV. To character- 
ize the spectra of these sources, we fitted each of the 7 bright- 
est X-ray sources with several different models (i.e., power-law 
(PL), thermal brem sstrahlung (TB), and blackbody (BB)) by 
using ICashI (Il979h statistics. Assuming all the X-ray sources 
within the half-mass radius are associated with the globular clus- 
ter M7 1 , we fixed the hydrogen column density at the value of 
A^H - 1.39 X 10^' cm - from optical extinction to attempt spec- 
tral fitting. 

In Table [H colum n 1 shows the Chandra source name given 
in lElsner et al.l (l2008l) . and column 2 lists the spectral model we 
used (we omitted the models that could not provide any physi- 
cally acceptable description of the observed spectra). Column 3 
gives the minimum number of counts used to group the spectral 
data for fitting, column 4 shows the best-fit photon index (F) or 
the temperature (keV), column 5-6 gives the C-statistic and the 
number of PHAQ bins, and the last column lists the unabsorbed 
X-ray flux in units of 10"'"* ergs s"' cm"^ in the energy bands 
0.3-8.0 and 0.5-2.5 keV. 

The X-ray spectra discussed here can help us to classify 
the faint X-ray sources. The brightest X-ray sources with Lx ^ 
10^^ ergs s"' in the energy band 0.5-2.5 keV and soft spectra 
(/0.5-2.0 kev//2.o-6.o keV ^ 1) are mostly qui escent low-mass X- 
ray binaries (qLMXBs JVerbunt et aLll2008 '). None of the X-ray 
sources in our sample shows this characteristic, and we con- 
clude that M71 does not contain this kind of binary systems. 
Cataclysmic variables (CVs) usually have hard spectra (power- 
law photon indices F < 2) and their X-ray luminosities are typ- 
ically between a few 10^" and a few lO"*^ ergs s"'. Most faint 

^ The number of source counts is a result of modelling the back- 
ground and the PSF, and represents a background subtracted expecta- 
tion value for the source counts within the detect cell. 

^ PHA: Pulse Height Analysis 



(Lx < 3 X 10^' ergs s ') sources with s oft spectra belong t o 
chromospherically active binaries (ABs; IVerbunt et al.l l2008l) . 
Looking at the brightest seven X-ray sources within the half- 
mass radius, we see that three have soft spectra (F > 2), three 
have hard spectra (F < 2), and one (s20) is borderline, with 
F ~ 2. These spectra suggest that s05, s08, and s29 might be 
CVs, AGN, or MSPs, while the softer spectra of s02, sl5, sl9 
and s20 may indicate that these are ABs. Definitive classifica- 
tions require optical identification, which we turn to now. 

3. Optical Observations 

Two fields located inside the half-mass radius of the globular 
cluster M71 were observed with the Wide Field and Planetary 
Camera 2 (WFPC2) on board the Hubble Space Telescope 
(HST) in 2000 and 2006. An image of the observations is shown 
in Figure [1] For these observations, the PC camera was cen- 
tered on the cluster center and the F336W (similar to U, hence 
t/336 hereafter), F439W (£439), and F555W (Vsss) filters were 
used. Exposure times were 800 sec in F336W (GO10524) and 
240 sec in F439W (G08118). Two exposure times correspond 
to the F555W filter are 80 sec and 63 sec for GO 10524 and 
G08118, respectively. GO10524 also contains F255W images 
which did not go deep enough to identify our targets. To esti- 
mate whethe r the n ondetections may be meaningful, we used 
iFerraro e7a n (12000) to estimate that CVs may be up to 3 mag- 
nitudes bright er (absolute magnitude) in F255W than V. Using 
ISeatonI (1 19791) . we estimate the extinction A255 = 7.0 for M71, 
and thus any CVs would be observed to be at least 3.2 magni- 
tudes fainter in F255 W than V. Using the WFPC2 exposure time 
calculator (ETC), we estimate that the brightest CV candidate in 
our WFPC2 field, s29, could attain a signal-to-noise ratio (SNR) 
of 1.4 in the F255W data if it showed the maximum F255W/y 
excess, which many of Ferraro et al.'s UV-selected objects do 
not. Therefore we did not discuss the F255W data further in 
this paper. The 5-cr limiting magnitudes of t/336, ^439, and 
for point sources are 21.09, 20.87, and 21.87, respectively. M71 
was also observed with the HST Advanced Camera for Surveys 
(ACS). The observations (GO10775) consist of F606W (Vsoe) 
and F814W (/gu) images covering the entire half-mass radius of 
the cluster (see Fig. [TJ. The exposure times for the F606W and 
F814W filters were 304 and 324 sec with 5-cr limiting magni- 
tudes of Vfjoe = 25.17 and /g^ = 23.91 for point sources. We 
note that the median value of the point sources with the signal to 
noise ratio S/N ~ 5 is used to define the 5-sigma limiting magni- 
tude. 

This section outlines the data reduction, photometry, and as- 
trometry of the HSTAVFPC2 and ACS images. 

3.1. Data Reduction and Photometry 

The HST/WFPC2 data obtained from the ESO archive were 
processed through the WFPC2 Associations Science Products 
Pipelin^ For each filter, single exposures were calibrated, in- 
cluding full bias subtraction and flat-fielding, and combined to- 
gether in order to remove the cosmic -ray events and coiTect 
the geometrical distortions 0. We also downloaded the archival 
HST/ACS drizzled images. Those images were combined from 
two Wide Field Channel ( WFC) images an d calibrated with 
MultiDrizzle package ( Koe kemoer et al.ll2002|) , which coiTected 
for geometric distortion and performed cosmic ray rejection. 

See http://archive.eso.org/archive/hst/wfpc2_asn/ 
^ see http://archive.stsci.edu/hst/wfpc2/pipeline.htrnl] 
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X-ray Source 


Model 


Grouping^ 


F/kT 


C-statistic/bins 


/o.3-8.o(,/b.5-2.5i^ 


s02 


PL 


8 


3.08+0.3 


12.11/9 


1.49(0.74) 




TB 


8 


0.62±0.1 


10.99/9 


1.26(0.82) 


s05 


PL 


30 


1.55±0.1 


8.20/9 


5.67(2.25) 




TB 


30 


9.96!|2 


10.33/9 


5.48(2.23) 




BB 


30 


0.69!«;^ 


59.72/9 


2.63(0.91) 


s08 


PL 


5 


1.70±0.5 


2.56/6 


0.74(0.32) 




TB 


5 


0.53 


20.18/6 


0.67(0.45) 


sl5 


PL 


5 


2.57±0.3 


9.34/1 1 


1 io//~\ ^\ 

1.18(0.63) 




TB 


5 


1 31+"-^ 


8.23/11 


0.84(0.57) 


sl9 


PL 


5 


2.61±0.4 


1.19/9 


1.14(0.60) 




TB 


5 


1 23+°'' 


2.77/9 


0.80(0.55) 


s20 


PL 


8 


1.92±0.4 


9.52/9 


1.71(0.83) 




TB 


8 


1 90+^-^ 

-oj 


8.58/9 


1.40(0.89) 




BB 


8 




0.84/7 


0.83(0.72) 


s29 


PL 


5 


1.29±0.4 


4.94/8 


1.12(0.36) 



" Minimum number of counts used to group the spectral data for fitting in XSPEC. 

* Unabsorbed X-ray flux in units of lO"'"* ergs s"' cm"^; energy bands are 0.3-8.0 (0.5-2.5) keV. 



Although M71 is a globular cluster, its stellar surface den- 
sity is not as dense as that of a typical globular cluster Even in 
the central region of M71, the average distance between stars 
is around 2", which is about 10 times larger than the typical 
FWHMs of WFPC2 and ACS cameras. Therefore, a simple aper- 
ture photometry method with the aperture correction is applica- 
ble to our data. We tested the flux measurement using several 
different psf-fitting photometry methods and the simple aperture 
photometry. We found that aperture photometry method had bet- 
ter signal-to-noise ratio and less magnitude error. Therefore, we 
decided to use aperture photometry to measure the fluxes of our 
data. 

For the data taken with WFPC2, we basically follo wed the 
instruction of aperture photometry described in Holtzma n et al.l 
(Il995l) . To deal with the PSF variances within each chip and 
between chips, we further separated the images into 4 and 9 
equal-size regions for the PC and WF chips, respectively, and 
performed the aperture photometry with aperture correction for 
each separated regions individually. We used an aperture with 
the size of (X'S in radius to measure fluxes, for all the objects 
with 3 <j detection found using the IRAF daofind task. We note 
that only less than 1% objects in each chip with separation of 
< (y.'S to their neighbors, so that using an aperture with the size 
of (y.'S does not suffer from the PSF overlapping problem. The 
local sky values were measured using an inner sky annulus of 
4 arcsecs with a width of 2 arcsecs, and the aperture correction 
value was calculated using the averages of the differences be- 
tween the magnitudes measured using apertures with sizes of 
(X'S and 4" in radius for 4 to 5 isolated stars in each separated 
region. The aperture correction valu e is 0. 1 1 +0.02 mag, wh ich is 
consistent with the value shown in 'Holtzman et al.l (Il995l) . The 
final output magnitudes, in the VEGAMAG system, were cor- 
rected for the appropriate zeropoints based upon the sensitivity 
infor mation in each header and the charge transfer efficiency ef- 
fect (lDolphinll2000l) . 

For the data taken with ACS, we perfor med the aperture pho - 
tometry based on the method described in ISirianni et al.l (l2005l) . 
The method is very similar to what we did for WFPC2. We also 
separated the ACS drizzled images into 9 equal-size regions to 
deal with the PSF variances. We used an aperture size of (X'S 
in radius to measure the fluxes with sky annulus from 4 arcsecs 
to 6 arcsecs. The aperture correction value is 0.08+0.01, which 



agrees with the values shown in lSirianni et alj ( l2005h . However, 
several optical counterparts, e.g., s08, sl9, suffer from the PSF 
overlapping problem since the distances between them and their 
neighbors are < I'.'O. In order to measure accurate fluxes for 
these counterparts, we first subtracted their neighbors by using 
the PSF generated from isolated stars which are close to the 
counterparts, and then performed the aperture photometry on 
these sources. By doing this, we can minimize the photometric 
effect from the PSF wings of neighbors. 

Com paring wi t h the photometry of M71 kindly provided 
by .Anderson et al.l ( |2008| ) for the ACS images and reported by 
iPiotto et alJ (12002) for the WFPC2 B-band and V-band images, 
their results are consistent with what we have obtained by us- 
ing aperture photometry, but have the main sequence and the 
giant branch with less noise. We therefore used their photometry 
in our study. In addition, for those possible optical counterparts 
undetected in their photometry, we used our own results for the 
magnitudes which have been corrected for the appropriate zero- 
point. 

The most informative of these diagrams are shown in Fig.|2] 
on which all stars located within the 95% confidence error cir- 
cles (see Section 2 and Table 1 of Eisner et al. 2008) of the 
Chandra source positions are indicated by red squares. Numbers 
have been assigned to all can didate counterparts corresponding 
to the 's' designation given in lEIsner et alj (l2008l) . with 'a', 'b', 
or 'c' appended if multiple potential optical counterparts exist. 

3.2. Astrometry 

To search for optical counterparts to the Chandra X-ray sources 
we aim to place both the X-ray and the optical frames onto the 
International Celestial Reference System (ICRS). We use this 
approach to improve the absolute pointing accuracy of Chandra 
and HST, a'6 an d I'.'O (Icr) respectively dAIdcroftet al.ll2000l: 
He ver et al.ll2004l) . 

Fo r the X-ray sources, the positions listed in lEIsner et"al] 
(l2008h are already on the ICRS. In this paper, we aim to tie 
the HST pointing to the ICRS by finding matches between 
stars appearing on HST images and stars with accurate posi- 
tions in th e Two Micron All Sk y Survey (2MASS) Point Source 
Catalog (iSkrutskie et al.ll2006h . On the basis of the HST point- 
ing information contained in each image header, we used the 
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Fig. 1. Chandra ACIS-S3 image of the globular cluster M71 within the energy range of 0.3 - 8.0 keV. The large circles are centered 
on the nominal center of the cluster and have radii Tc (inner circle), (middle circle), and Ir^ (outer circles). The small circles 
show the positions of the 63 X-ray sources within vmii < 2rh jElsner et al ]l2008h . Straight lines mark the nominal boundaries of the 
ACIS-S3 and ACIS-S2 CCDs, with most of the figure falling on S3 and the lower left hand portion on S2. The field of view of the 
HST ACS (GO10775) marked by the black square covers the entire half-mass radius. The green (dashed) and red (solid) polygons 
are the field of view of the HST WFPC2 for GO 10524 and G08 118, respectively. 



WCSTools/Zmwcj^task on each corrected image to do the cross- 
correlation. The resulting positions were matched to those stars 
from the 2MASS catalog. There are hundreds of 2MASS stars 
within each HST image. By using those 2MASS stars as refer- 
ence, the astrometric solution yielded root-mean-square resid- 
uals of 0'.'057 in right ascension (RA) and (y.'065 in declination 
(Dec) relative to the 2MASS astrometry for the ACS images. The 
resultant solution gave the residual errors of 0'.'068 and 0'.'148 
in RA and a.' 066 and a.' 149 in Dec relative to the 2MASS as- 
trometry for the PC and WF images, respectively. The final un- 
certainties of the optical source position in RA and Dec are the 
root of the square sum of the uncertainty of the astrometry in 
2MASS and HST image alignment and the general uncertainties 
of 2MAS S point source astrome try of typically ~ CK.'l relative to 
the ICRS (ISkrutskie et al.ll2006h . 



4. Source Identification and Classification 

To obtain optical identifications for the X-ray sources we use the 
precise astrometry described in Section 3. We search for optical 
counterparts within the 95% Chandra error circle of the source 
positions (see Table 1 of lElsner et alj|2068l) . which includes the 
positional uncertainty of X-ray sources reported by the wavelet 



^ See |http://tdc-www.harvard.edu/software/wcstools/index.html| 



source detection algorithm, the uncertainty in the X-ray bore- 
sight correction, and the uncertainty in the optical astrometry. 
Within the half-mass radius of M71, there are 29 sources de- 
tected by Chandra and we suggest optical counterparts based on 
positional coincidence alone to 25 of them. In case of multiple 
sources inside the X-ray error circle, we include all the candi- 
dates. The results of each candidate optical counterpart are sum- 
marized in Table |2] and finding charts are shown in Fig. [3] and 

Fig. a 

The first step in classifying faint (Lx ^ lO^'^''' ergs s"') X- 
ray sources is to study their X-ray properties, e.g. their X-ray 
luminosity and spectral behavior (see §2). The second step in 
the identification process can be made when the sources have 
other information coming from different wavelengths. The coin- 
cidence between accurate radio timing positions of millisecond 
pulsars and the positions of X-ray sources can provide reliable 
identification. In the optical band, the color-magnitude diagrams 
(CMDs) of globular clusters have been studied for a long time 
because they reflect the fundamental properties of these stars and 
the evolutionary stage of the globular clusters. We extract further 
information from the locations of the optical stars in the CMDs 
of Fig. |2] CVs usually lie much bluer than the main sequence 
stars in the (V, U-V) and (V, B-V) CMDs while the X-ray ABs 
may be located on or slightly above the main sequence or on 
the giant branch. Stars on the main sequence in the V vs. V-I 
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Fig. 2. Color magnitude diagrams (CMDs) for all the sources detected in the WFPC2 and ACS field of view. The HST candidate 
counterparts matched to the X-ray sources are indicated by red squares. We note here sl4b in the (V, V-I) CMD is plotted as a 
leftward-pointing arrow since its /814-band magnitude is far below the 5-sigma hmiting magnitude and with a large magnitude error. 



CMD cannot be clearly classified: they could either be CVs or 
ABs since their optical flux is dominated by the donor stars or 
brighter stars of the binary systems. 

The ratio of X-ray to optical flux is also useful to distin- 
guish CVs from X-ray ABs (see iBassa e"tan f2004). In Fig. |5] 
we show the X-ray luminosity as a function of the abso- 
lute magnitude for low-luminosity X-ray sources from 47 Tuc, 
NGC 6397. NGC6752, M4 , NGC288. M55 NGC 6366 (data 
from iGrindlav et al.' '2001a; Edmonds et al. 2003a; 'Cool et al.' 
1998; Grindlayetal. JOO lb; Taylor etal. 20()^1; Pooley et al 
2OO2I: iBassa et al.1 i2004 iKong et alJ l2006i iBassa et al.1 feOOSh . 
and M71. The large symbols in this figure indicate the X-ray 
sources with possible optical counterparts in the field of view 
of the Chandra observation of M71, while the smaller symbols 
show classified objects found in other clusters. We note that the 
absolute magnitudes and X-ray luminosities for the sources in 
the observations are computed under t he assumption that they 
are cluster members. As discussed in lElsner et a n (I2008h . we 
caution that ~ 40% of the 29 X-ray sources within the half-mass 
radius are background or foreground objects. 

Now we turn to those stars unrelated to the globular clus- 
ters. Foreground stars are likely to have counterparts not on 
the main sequence, have soft spectra, and have low fx/fo ra- 
tios (lKrautteretalJll999[) . For background active galactic nu- 
clei (AGN), they will also have counterparts not on the main 



sequence, but with hard spectra; their fx/fo ratios can be high 
(Krautteret al. 1999). However, the AGN won't necessarily be 
detected at all; in some cases, the only object in an error circle 
may be a cluster main-sequence star that is not related to the 
X-ray source. 

We first consider those X-ray sources with only one sug- 
gested counterpart in the Chandra error circle. S08 is a known 
miflisecond pulsar, PSR J1953-I-1846A = M71A, with a spin 
period of 4.89 ms. It is in a 4.24 hr eclipsing binary sys- 
tem with a low-mass (> 0.032 Mq) companion. It was dis 



tem witn a low-mass (> U.UjZ Mq) companion, it was ais- 
covered with Arecibo dRansom et alj|2003l 120051; iHessels et al.l 



I2007h and its radio emission is partially eclipsed in the orbital 
pha se interval .18 - 0.36 for approximately 20% of each or- 
bit (Hessels et al. 2007). The X-ray counterpart was detected by 
Eisner et al. (2008 ). Within the Chandra error ckcle we find a 
possible optical counterpart of this pulsar in the Veoe and hi4 
band. The candidate optical counterpart to s08 lies on the main 
sequence in our V-I CMD with an absolute magnitude My ~ 8.5, 
which implies that it has a mass of about 0. 5 Mq. The radio tim- 
ing indicates a minimum mass of 0.03 Mo (IHessels et alj|2007h 
for the companion star In order to allow for such a massive 
companion (i.e. ~0.5 Mq), conceivably the orbit could be seen 
nearly face-on (within 4 degrees). However, in that (extremely 
unlikely) case we do not expect regular radio eclipses, as are ob- 
served. M71 A's radio properties are very similar to those of other 
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Fig. 3. 5 " X 5 " finding charts for candidate optical counterparts within the half-mass radius of M7 1 . The finding charts are obtained 
from the HST/ ACS 7606 images. The 95% confidence uncertainties on the Chandra positions are overlaid on these charts, while 
the candidate counterparts are indicated with an arrow. The finding charts are set with their dynamic range. The greyscale of these 
images is chosen to enhance the visibility of the candidate counterparts. All images have North to the top and East to the left. 
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Fig. 4. 5" X 5" finding charts for candidate optical counterparts outside the half-mass radius of M7 1 . The finding charts are obtained 
from the HST/ACS Va)6 images. The 95% confidence uncertainties on the Chandra positions are overlaid on these charts, while 
the candidate counterparts are indicated with an arrow. The finding charts are set with their dynamic range. The greyscale of these 
images is chosen to enhance the visibility of the candidate counterparts. All images have North to the top and East to the left. 



very low-ma ss binary pulsars such as PSR J1701-3006E (M62E, 
lFreirel2005h . and therefore we conclude that M71 A's companion 
is not this star. Although this star's position agrees within (X'l 
with M71A's position from radio timing (I. Stairs 2009, private 
comm. ), transferred with O*.' ! accuracy (Icr) onto the 2MASS 
frame dSkrutskie et alj|2b06l) , this could be coincidence due to 
the crowding in this field (Fig.[3]i; alternatively, M71A could be 
a hierarchical triple system. Future radio timing may determine 
this. 

The star in the error circle of s02 is nearly located on the 
main-sequence turn-off point (MSTO) in the (V, V-I) CMD of 
Fig. |2] and slightly below the subgiant branch, which is simi- 
lar to the "red stragg ler" active binar i es seen as X-ray sourc es 
in other clusters ( Albrow et al J 120011: lEdmonds et"al]|2003allbl: 
[Bassa et al. 2008). Its ratio of the X-ray to optical flux locates 
in the region of ABs in Fig.|5] Since s02 has a soft s pectrum and 
shows significant time variability in the X-ray band CElsner et alj 
120081) we suggest that s02 is a chromospherically AB and its 
temporal variation can be explained as flaring on the coronally 
active star. S04 and sl8 are also believed to be in the same group 
as ABs since both of them are located slightly above the main- 
sequence turn-off point, have soft X-ray spectra, and have lower 
X-ray to optical flux ratios. 

SOS is the brightest X-ray source within the half-mass ra- 
dius of M71. It is worth noting that s05 is the only X-ray source 
detected with ROSAT jPanzera et al.1 120031: [Eisner et al.ll2008l) 
inside the half-mass radius. Its optical counterpart is bluer than 
the main-sequence, and it has a relatively high X-ray luminos- 
ity (Lx ~ 4 X 10^' ergs s"^). It is unlikely to be an AB. Its 
X-ray spectrum is too hard to consider it as a quiescent low- 
mass X-ray binary (qLMXB). SOS gives a bremsstrahlung tem- 
perature consisten t with ~ 10 keV, as appropriate f or luminous 
magnetic CVs (E racleous et alj|l99ll: lMukaill2003b . During the 
S2.4 ks observation time, it is not consistent with being steady at 
99.9% confidence. We suggest that a CV interpretation is plau- 
sible. In addition, sOS has a high value of log(/x//o) ~ 0.37 
and the blue color, which implies this source could be a back- 
ground AGN. However, the power-law fit of its X-ray spectrum, 
with photon index r = l.SS + 0.1, might be considered as arising 
from the intra-binary shock formed due to interaction between 
the relativistic pulsar wind and material from its companion star. 
An irradiate d main-sequence com panion could be this blue; e.g. 
47Tuc W tBogdanov et alj l200Si) . MSPs with main-sequence 



companions of the mass of ~ O.SMq have not yet been detected, 
but may well be hidden from r adio detection by clouds of ionized 
gas from the companion (e.g. lFreire et al.l l2004). Therefore, we 
cannot rule out the interpretation that it is a binary MSP system 
though this unusual scenario must be judged unlikely. 

The candidate cluster counterpart to s29 has ultraviolet ex- 
cess with respect to the main sequence (Fig. |2]l and has a high 
X-ray to optical flux ratio. The source can be well fitted with a 
power-law model with a photon index of F = 1 .29 + 0.4, and 
its X-ray luminosity is ix,o.3-8.o keV ~ 2.1 x 10"" ergs s"'. Its 
U-V color is far too blue to be an AB while the optical color is 
redder, almost on the main sequence. That indicates s29 is a CV 
with two spectral components, a blue disk and a red companion 
star. 

Source s IS is a good AB candidate since there is no evidence 
for a blue color in the VI CMD and it has a soft X-ray spectrum 
(/0.5-2.0 kev/ /2.0-6.0 keV ^ !)■ On the other hand, without the in- 
formation from the U-V or B-V color, a CV interpretation is 
still plausible. Its relatively high value of log(/x//o) ~ -0.79 
suggests that slS could be a CV, although it doesn't rule out an 
AB. 

The star in the error circle of s20 lies on the giant branch 
and has a soft X-ray spectrum, which gives strong evidence that 
it is a chromospherically active binary containing a giant star 
i.e., a RS CVn sys tem). Its temporal variation in the X-ray band 
Eisner et al can be explained by magnetic activity. Since 

s07 and sl2 are located on the giant branch and have relatively 
low X-ray to optical flux ratios, we believe that they are likely 
RS CVn systems as well. 

The optical counterparts associated with those X-ray sources 
having lower photon statistics, sOl, s06, and s23, are located on 
the main-sequence and have lower X-ray to optical flux ratios. 
We then consider that all of them may be X-ray ABs. S22 ex- 
hibits rather interesting colors. In the V-I CMD it is on the main 
sequence, but as we shift to progressively bluer colors its color 
gets redder and redder while in the U-V CMD it is way off the 
main sequence. Thus we suggest that it is either a foreground or 
background source, not associated with M71. 

We turn now to the sources with more than one possible 
counterpart in the error circle. We find two or three possible op- 
tical counterparts within each of the Chandra error circles for 
s03, sl4, sl9, s27, and s28. S03a and s03b both fall on the main- 
sequence in the (V, U-V), (V, B-V), and (V, V-I) CMDs. Their 
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Table 2. Optical counterparts to Chandra X-ray sources within the HST/ ACS field of view 



X-ray 


Offset From CXO position^ 


t/336 


S439 


V555 


^606 


^814 


log(Lxt 


log(/x//oL 


FAP^ 


Classification^ 


Source 


ARA (") 


ADec (") 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


ergs s"' 








Inside the half-mass radius 


sOl 


0.22 


0.49 


. . 






20.85(0) 


19.94(2) 


29.85 


-2.02 


26.06 


AB? 


s02 


0.0 


-0.05 








17.38(1) 


16.32(1) 


31.18 


-2.08 


1.08 


AB 


s03a 


0.18 


-0.69 


18.68(6) 


18.35(4) 


17.78(4) 


18.33(1) 


17.56(0) 


29.95 


-2.93 


7.01 


AB? 


s03b 


-0.50 


-0.29 


19.32(9) 


19.09(4) 


18.44(4) 


18.90(1) 


18.11(0) 


29.95 


-2.70 


10.36 


AB? 


s03c 


-0.57 


0.37 




X 


X 


23.8(1) 


22.81(6) 


29.95 


-0.76 


37.35 


CV/AGN? 


s04 


0.07 


-0.02 








16.92(1) 


16.01(1) 


30.63 


-2.82 


0.73 


AB? 


s05 


0.06 


-0.05 








22.42(5) 


21.79(0) 


31.61 


0.37 


11.98 


CV/AGN? 


s06 


0.04 


-0.08 








17.74(1) 


16.91(1) 


30.15 


-2.96 


2.67 


AB? 


s07 


0.33 


-0.11 


16.26(2) 


15.16(5) 


14.30(5) 


14.86(0) 


13.89(0) 


30.31 


-3.96 


0.25 


AB? 


s08 


0.24 


-0.17 


X 


X 


X 


21.48(2) 


19.96(1) 


30.78 


-0.84 


11.31 


MSP 


slO 


0.27 


0.05 


X 


X 


X 


25.25(6) 


23.4(1) 


30.37 


0.26 


21.96 


CV/AGN? 


sl2 


0.28 


-0.10 


16.40(2) 


15.60(3) 


14.93(3) 


15.55(0) 


14.82(0) 


29.88 


-4.11 


0.72 


AB? 


sl4a 


-0.29 


0.54 


X 


X 


X 


23.47(6) 


20.87(2) 


30.03 


-0.80 


33.60 


F? 


sl4b 


-0.04 


-0.18 


X 


X 


X 


25.8(1) 


28(2) 


30.03 


0.12 


48.35 


CV/AGN? 


sl4c 


0.18 


0.36 


X 


X 


X 


26.4(2) 


23.46(6) 


30.03 


0.36 


51.13 


CV/AGN? 


sl5 


0.10 


0.06 








21.06(6) 


19.98(4) 


31.0 


-0.79 


10.37 


CV/AB? 


sl6 


0.15 


-0.09 


19.9(1) 


18.62(4) 


17.55(3) 


17.95(1) 


16.88(0) 


30.15 


-2.88 


3.44 


F? 


sl8 


0.08 


-0.08 


17.25(3) 


16.96(6) 


16.28(6) 


16.69(1) 


15.92(1) 


30.77 


-2.77 


0.58 


AB 


sl9a 


0.08 


-0.12 


18.20(5) 


17.79(3) 


16.99(3) 


17.39(0) 


16.50(0) 


30.95 


-2.30 


1.16 


AB? 


sl9b 


-0.28 


0.10 


19.52(9) 


18.93(4) 


18.32(4) 


18.83(1) 


18.05(1) 


30.95 


-1.73 


4.62 


AB? 


sl9c 


0.18 


0.36 


X 


X 


X 


21.48(1) 


20.04(2) 


30.95 


-0.66 


11.33 


AB? 


s20 


0.08 


-0.07 


16.16(2) 


15.21(3) 


14.18(4) 


14.65(0) 


13.61(0) 


31.18 


-3.17 


0.16 


AB 


s21 


0.17 


-0.15 


X 


X 


X 


27.09(9) 


26.5(3) 


29.99 


0.61 


42.03 


CV/AGN? 


s22 


0.14 


-0.15 


19.01(7) 


18.22(5) 


17.52(3) 


18.03(0) 


17.22(1) 


30.02 


-2.98 


5.48 


F? 


s23 


0.32 


0.09 


18.11(5) 


17.75(3) 


17.19(3) 


17.73(1) 


16.73(1) 


30.22 


-2.91 


2.49 


AB? 


s24 


0.21 


0.05 


X 


20.3(1) 


19.13(4) 


19.4(1) 


18.16(1) 


30.02 


-2.43 


9.81 


F 


s25 


0.08 


0.06 


X 


X 


X 


25.49(8) 


23.3(1) 


30.03 


0.01 


30.34 


CV/AGN? 


s26 


-0.07 


-0.08 




19.3(2) 


18.33(7) 


18.63(7) 


17.50(3) 


30.66 


-2.10 


4.88 


F 


s27a 


0.19 


0.17 


16.06(2) 


15.15(4) 


14.23(4) 


14.68(0) 


13.70(0) 


30.22 


-4.12 


0.24 


AB? 


s27b 


-0.57 


0.06 


19.18(8) 


18.41(3) 


17.73(3) 


18.19(1) 


17.45(1) 


30.22 


-2.72 


4.09 


AB? 


s28a 


0.35 


-0.06 


X 






20.67(0) 


18.81(0) 


29.99 


-1.96 


17.96 


F? 


s28b 


-0.15 


0.44 


X 






20.28(1) 


19.28(1) 


29.99 


-2.11 


16.12 


AB? 


s29 


0.05 


0.06 


21.0(2) 




20.9(2) 


21.44(1) 


20.76(1) 


30.85 


-0.79 


11.20 


CV 


Outside the half-mass radius 


s41 


0.13 


-0.12 








21.53(1) 


20.77(1) 


30.69 


-0.91 


14.47 


CV/AGN? 


s49a 


-0.24 


0.09 




19.87(7) 


18.49(3) 


18.74(1) 


16.72(1) 


30.13 


-2.59 


9.43 


F 


s49b 


-0.31 


-0.38 




X 


X 


20.93(2) 


18.35(1) 


30.12 


-1.71 


21.64 


F 


s54a 


0.68 


0.20 








15.36(0) 


14.46(0) 


30.62 


-3.45 


0.65^ 


AB? 


s54b 


-0.23 


-0.05 








21.46(1) 


20.72(1) 


30.62 


-1.01 


19.76'' 


CV/AGN? 



" The position offsets of sources detected in the V^oe band are given relative to their Chandra positions (see Eisner et al. 2008 for details) 

* Lx'. unabsorbed X-ray luminosity in the energy band is 0.5-2.5 keV (see Table 5 of Eisner et al. 2008^) 
Ratio of X-ray to optical (Keoe) flux is computed by using log(/x//o) = log/x -1- 0.4^606 + 5.07 ; fx is derived in the 0.5-2.5 keV band. 
FAP: False Alarm Probability (%). The probability of a chance coincidence in % is calculated f rom Pcoinc = »^)^ ?'q^ /^search (see §6.1 
of lElsner et al.ll2008l) . The search area, Asearch, is the central region of M71 with the radius of 2' (see [Anderson et al]|2008l) and the considered 
magnitude, m, is Vjoe band magnitude. 

'' MSP: millisecond pulsar; CV: cataclysmic variable; AB: X-ray active binary; AGN: active galactic nucleus; F: foreground source; ?: all 
cases where we are not certain about the nature of the X-ray source. 

' Sources lying outside the field-of-view of the WFPC2 chips do not have F336W, F439W or F555W data 

" Under detection limit. 

Since the source is outside the searching region of the 2'-radius central circle dAnderson etai]|2008h , the Aseaich for the probability of a 
chance coincidence is calculated with the HST/ACS field of view, i.e. 202" x 202". 



colors and low X-ray to optical flux ratios suggest that either 
s03a or s03b is a chromospherically AB. However, the blue color 
and the relatively high X-ray to optical flux ratio of s03c indi- 
cates that it is a CV or a background AGN. S14a and s28a are 
located far from the main sequence, suggesting that they do not 
belong to M71, while s28b located on the main sequence could 
be an active binary system since it has the relatively low X-ray 
to optical flux ratio. In Fig. [3] we find two additional possible 



optical counterparts, sl4b and sl4c, within the Chandra eiTor 
circle of sl4, which are fainter than the 5-cr limiting magnitudes 
of Veoe and /814. Based on the blue color and the high value of 
fx/fo, si 4b could be either a CV or a background AGN. The 
optical-faint source sl4c is located near the downward-extended 
part of the main sequence, suggesting that it could be a main- 
sequence star. However, it is located in the Lx vs. My diagram 
in a region where no authentic cluster members have been found 
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if we compute its X-ray luminosity and absolute magnitude un- 
der the assumption that it belongs to M7 1 . We then rule out the 
AB interpretation. Therefore, due to its high X-ray to optical flux 
ratio of this source, sl4c is considered as either a good candidate 
for background AGN with a optically faint object inside the er- 
ror circle not related to the X-ray source or a CV candidate with 
a secondary star that dominates the optical flux in the (V, V-I) 
CMD. There are three optical counterparts within the Chandra 
error circle of sl9. The position of sl9a is near the MSTO point 
and slightly below the subgiant branch, which is similar to the 
case of s02, while sl9b is located on the main sequence. Both of 
their X-ray to optical flux ratios are located in the region that is 
primarily populated by ABs (Fig.|5]l. For the third optical coun- 
terpart, si 9c, its location in the VI CMD and soft X-ray spectrum 
indicate that si 9c is an active binary as well. However, without 
the information from the UV color we cannot eliminate the CV 
interpretation due to its relatively high X-ray to optical flux ra- 
tio, which is similar to the case of sl5. We suggest that sl9a is 
the most likely counterpart, as re d stragglers are very often asso- 
ciated with X-ray sources (e.g. iHeinke et alj|2005l) . According 
to the positions of two possible optical counterparts to s27 in the 
CMDs and in the absolute magnitude vs. X-ray luminosity dia- 
gram, we believe that either s27a or s27b is likely to be an active 
binary. 

The optical counterparts of si 6, s24, and s26 are located fur- 
ther above or to the right of the main sequence than the binary 
sequence. Hence we believe that they are foreground objects and 
unrelated to M71. The positions of two faint optical counter- 
parts to slO and s25 in the (V, V-I) CMD and their relatively 
high X-ray to optical flux ratios are very similar to the case of 
sl4c so that slO and s25 could be either CVs or background 
AGN. The highest X-ray to optical flux ratio among 39 possible 
counterparts and blue color suggests that s21 is the most likely 
background AGN although we cannot eliminate the interpreta- 
tion of a CV. Furthermore, the regions of the X-ray sources, s09, 
si 1, sl3, and sl7, were also observed with the HST/WFPC2 and 
ACS, but we do not find any optical counterparts inside their 
Chandra error circles. If we set the S-cr limiting magnitude of 
V(,Q6 = 25.17 as the upper limit for these sources, their X-ray to 
optical flux ratios fall on a range of ~ 0.2-0.6, which are higher 
than the highest log( /x//o) value known fo r an AB in a cluster, 
e.g. W64 in 47Tuc jEdmonds et al.ll2003al) . Therefore, an AB 
interpretation can be rejected. If we take their X-ray colors into 
account, sll and sl7 are located near the bottom-right corner 
and close to the po sition of s40 i n the X-ray color-color diagram 
shown in Fig. 4 of lElsner et al.l (l2008l) . which implies they have 
very hard spectra with over half of their counts above 2 ke V. This 
infers a high intrinsic A^h, which strongly suggests that these are 
background AGN. For the other three sources with medium X- 
ray colors, we then tentatively classify them all as CVs, MSPs, 
or background AGN, though AGN are probably the most likely 
category. 

Outside the half-mass radius of M71, we find 6 optical coun- 
terparts to X-ray sources, s41, s42, s49 and s54 in the ACS field- 
of-view. The candidate counterpart to s42, located on the edge of 
the ACS, is saturated in the optical band, which prevents us from 
obtaining a reliable magnitude of this optical source or search- 
ing for any other faint optical sources inside the Chandra error 
circle. There are 2 possible optical counterparts to s49. Both of 
them are located far from the main sequence, hence we suggest 
that they are not associated with M7 1 . Inside the error circle of 
s54, the brighter object s54a is on the giant branch while the 
fainter one, s54b, is located blueward of the main sequence in 
the CMD and has a relatively high X-ray to optical flux ratio. 
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Fig. 5. X-ray luminosity as a function of absolute magnitude, for 
low-luminosity X-ray sources in globular clusters. Five types of 
X-ray sources are shown, qLMXBs (diamonds), MSPs (crosses), 
CVs (triangles), ABs (stars), and unclassified sources (squares). 
The larger and numbered symbols in this figure correspond 
to the optically identified X-ray sources in the field of view 
of the Chandra observation of M71, where we compute abso- 
lute magnitude and X-ray luminosity under the assumption that 
the sources are associated with M71. The smaller symbols in 
this figure indicate objects found in other clusters, i.e. 47 Tuc, 
NGC 6397, NGC 6752, M4, NGC 288, M55, and NGC 6366. We 
note that ambiguous sources coming from other clusters were 
discarded in this figure. The dashed line of constant X-ray to op- 
tical flux ratio given b y log ix,o.5-2.5 kev(ergs s"') = 34.0 - 0.4 
My (after iBassa et al]i2004]) roughly separates CVs from ABs. 

suggesting that s54b might be a CV. However, s54b lies on the 
spikes produced by s54a, which prevents us from obtaining an 
accurate magnitude for s54b. S4rs color is bluer than the main 
sequence, and its X-ray to optical flux ratio is higher than that of 
an AB. We then suggest that s41 is a CV candidate. Furthermore, 
a background AGN scenario is also plausible for s41 and s54b 
due to their blue colors, high log(/x//o) values, and their loca- 
tions outside the half-mass radius of M71, where they are more 
likely to be background sources. 

5. Summary and Discussion 

In summary, we find one certain CV (s29), seven possible candi- 
date CVs (s05, slO, sl4, s21, s25, s41, and s54), and two certain 
ABs (s02 and s20) and 12 good candidate ABs (sOl, s03, s04, 
s06, s07, sl2, sl5, sl8, sl9, s23, s27, and s28) in the globular 
cluster M7 1 . Some of our candidate CVs (and/or candidate ABs) 
might be MSPs in binary systems or AGN, which often (but not 
always) show blue colors. 

To interpret our results, understanding how many of our ob- 
jects are likely false matches will be critical. We calculated the 
expected number of false matches in several ways. First, we 
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shifted all our X-ray source positions by 9" and 18" (somewhat 
arbitrary, but chosen to be larger than the largest uncertainties) in 
four directions, and searched for matches against the Veoe frame. 
From this exercise, we expect 1 l^g false matches among our 39 
possible matches, indicating that > 70% of the 34 total X-ray 
sources in our field of view have a true match. By chance, then, 
70% of our false matches should occur with sources which have 
a true match — suggesting that ~ 7.7 sources should have two 
possible optical counterparts. We see seven sources that have two 
or more possible Vgofe counterparts, which is nicely consistent. 

In order to know which sources are more likely to have 
false matches, we calculated the probability of a chance coin- 
cide nce in % shown in Col. 10 of Table |2]by using Eq. (3) 
from lElsner et all ( |2008|) (see also Verbunt et al. 2008). Within 
the half-mass radius, there are 6 sources with a false alarm 
probability smaller than 1%, which indicates that those associ- 
ations between X-ray sources and optical counterparts have a 
>99% confidence level. Adding up the false alarm probabilities 
gives a total expected number of false matches of ~ 5, which 
is consistent with the the expectation of 11;^^ false matches 
above. Among the 10 X-ray sources with optical counterparts 
and Lx,o.5-6.o kev > 4 X 10"^" ergs s"', adding the false alarm 
probabilities indicates that roughly one of them is expected to 
be a false match. (Note that we believe s08 to be a false match, 
but to be a true cluster member; and that sl9 has three potential 
counterparts.) 

Inside the half-mass radius of M71, we find 14 X-ray sources 
with Lx,o.5-6.o kev > 4 X lO-'" ergs s"', of which 10 have optical 
counterparts. Assuming that all X-ray sources outside the half- 
mass radius are fore- or background sources, we can estimate 
that ~ 3.7^j g X-ray sources among 14 are unr elated to M71 . 
The error quote here is from the Poisson statistic (lGehrelsll986l) . 
This is consistent with our estimate above that we have identified 
true optical counterparts for ~ 9 X-ray sources. 

iPoolev et alj (l2003b have shown that for 12 globular clus- 
ters observed by Chandra, the number of globular cluster X-ray 
sources which are above the lower limit of 4x 10-"' ergs s"' (0.5- 
6keV) can be approximately linearly fitted with the predicted 
stellar encounter rate P oc Po^^^rc^ = F, where F is referred to as 
the collision number (I Verbunlll2003h . Here po is the central den- 
sity of the cluster, and is the core radius. In order to examine 
if M71 fits this relation we compare its number of X-ray sources 
and its collision number F with those of some other clusters, 
NGC6266, 47Tuc, M28, M4, NGC6366, M55, and NGC288 
(see Fig. |6ll, using the parameters listed in Table [3] We note 
that the core-collapsed globular clusters are not considered in 
our study since their core parameters are generally uncertain, in- 
troducing strong uncertainties into interaction rates derived from 
those parameters. The encounter number for M71 is ~ 230 and 
~ 10 times smaller than those of 47 Tuc and M4, respectively. 
iPooley et al.l (l2003h reports 41+2 sources a bove the lower lumi- 
nosity limit in 47 Tuc, which are revised by lHeinke et al.l (l2005h 
to 63±4 (for a distance of 4.5 kpc); the uncertainty is due to 
the estimated number of background sources. Thus, if the num- 
ber of sources scales with the encounter rate, the presence of 
~ 10 + 3 sources with Lx,o.5-6.o keV > 4 x lO"*" ergs s"' in M71 
is a very significant overabundance, even if we take into account 
the errors due to Poissonian fluctuations. The same conclusion 
is reached on the basis of comparison with any other globular 
clusters listed in Table[3]except for M55 and NGC 288, in which 
the number of the X-ray sources is also in excess of the predicted 
value. This indicates that most of the sources in M71, M55, and 
NGC 288 are not formed via stellar encounters. 



47 Tuc 
5 



M,4 



Collision Number r 



Fig. 6. Number of globular cluster X-ray sources (N) with 
Lx,u.5-6.o keV ^ 4 X lO""' crgs s"' inside the half-mass radius vs. 
the normalized collision number F of the cluster. The dashed line 
indicates the best power-law fit of the data points from 47 Tuc, 
NGC 6266, M28, M4, and NGC 6366 (cf. Table©. The relation 
between N a nd F, i.e. N oc p"-74±o .03^ jg consistent with the result 
derived from lPoolev et al.l (l2003h . 



As suggested bv' Verbunt* (l2002h . ABs are most likely primor- 
dial binaries, and thu s to first order thei r numbers should scale 
with mass. Following iKong et al.1 (l2006l) . we calculated the half 
masses with 10""-^'*^'', assuming the visual mass-to-light ratio is 
the same for all clusters listed in Table |3] M71 has the lowest 
half-mass, containing only ~ 30% of the mass within the half- 
mass radius of M 4. Scaled by mass, the predicted number of 
ABs with Lx.{).5-6kev > 4 X lO-^" ergs s"' in M71 should be 
similar to that of NGC 6366 and smaller than those in any other 
cluster shown in Table |3] but this is in contrast to our results. 

The scahng of source number with the collision number for 
the sources with Lx,o.5-6 keV > 4x10-"' ergs s"' suggests that CVs 
are mostly made via stellar encounters (Poolev et al. 2003). If we 
assumed that all of the CVs in M71 were formed dynamically, 
we would not expect to find more than one CV by scaling with 
the encounter numbers from any other cluster liste d in Table [3l 
However, most of the globular clusters studied by Poolev et al.l 
(2003) have high F numbers and many dynamically produced 
CVs. It is reasonable to suspect that primordial CVs may domi- 
na te in the low-d ensity clusters. According to the computations 
bv lDaviesI (Il997h . a cluster core with a star density of 1000 pc"^ 
allows most of the CV progenitors to evolve into a CV. This 
could explain the existence of at least one (and several candi- 
date) CVs within the half-mass radius of M71. 

It is interesting to mention that in M7 1 ~7 optical counter- 
parts to Chandra X-ray sources are classified as potential RS 
CVn systems, in which X-ray emission is produced primarily 
in (sub)giant flare outbursts. A couple possible RS CVn X -ray 
sources hav e been identified in 47 Tuc dHeinke et al.ll2005h and 
Lj Centauri ( Cool et al.ll2002l) . but, beside s M71, on ly the low- 
density clusters M55 and NGC 6366 (Ba ssa et al] 1^2008) have 
significant fractions of X-ray sources identified as possible RS 
CVns. Considering the even lower density case, the total X-ray 
luminosity o f the old open cluster M 67 is dominated by binaries 
with giants ( van den Berg et all 12004). From the ROSAT cen- 
sus, the X-ray emission of most globular cluster s per unit mass 
is lower than that of the old open cluster M67 (IVerbuntil200ll 
120021) . There are three possible explanations for this: a) M67 is 
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(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Cluster 


logpo 


rc 
(") 


d 
(kpc) 


My 


r 


Mh 


Source 


Background 


Member 


IDs. 


Reference" 


NGC 6266 


5.14 


10.8 


6.9 


-9.19 


37.07 


8.24 


51 


2-3 


48-49 


2' 


1,2,3 


47 Tuc 


4.81 


24.0 


4.5 


-9.4 


24.91 


10.0 


79 


-16 


63±4 


53-63 


4 


M28 


4.75 


14.4 


5.6 


-8.33 


11.29 


3.73 


26 


2-3 


23-24 


T 


1,5,6 


M4 


4.01 


49.8 


1.73 


-6.9 


1.0 


1.0 


6 


1-3 


3-5 


5 


1,7 


M71 


3.05 


37.8 


4.0 


-5.6 


0.11 


0.30 


14 


1-7 


10±3 


4-9 


8 


NGC 6366 


2.42 


109.8 


3.6 


-5.77 


0.08 


0.33 


5 


2-5 




~1 


9 


M55 


2.15 


169.8 


5.3 


-7.6 


0.18 


1.82 


16 


5-12 




2-4 


9 


NGC 288* 


1.80 


85.0 


8.4 


-6.7 


0.03 


0.83 


11 


4-11 


4- 


2-5 


10 



Note: Cols. 2-5 give the values for central density (po), core-radius (r^), distance (d) and absolute visual magnitude (My) originate from' HarrisI 
version of February 2003). For M4, the values of po and My are computed for the distance and reddeni ng of Richer et al . (1997). Cols. 
6 and 7 are the collision number which is computed from F = p^^ and the half-mass from Mh oc lO-°'^'^v dKong et alj|2006h . Values for F 
and Mh are normalized to the value of M4. Col. 8 shows the total number of sources detected within the half-mass radius. Col. 9 is the number 
of expected fore/background sources. Col. 10 gives the number of expected cluster members plus error. Col. 11 shows the number of X-ray 
sources (with Lx,o.5-6.o kcv > 4 x 10'" ergs s"') which have optical or/and radio counterparts associated with the cluster, or spectrally confirmed 
qLMXBs. The last column gives the reference paper. These globular cluste rs are ordered on the cen tral density. 

a. For each cluster , th e basic data for this table were extracted from: 1. Poolev et al. (2003); 2. Cocozza et al. (2008); 3. Trepl (2007); 4. 
iHeinkeetal .1 ( 120051) ; 5. lBecker et all ( 120031) ; 6. Becker & Hui (2007); 7. Bassa et al. (2004); 8. this work; 9. Bassa et al. (2008); 10. Kong et al] 
( I2006h . 

b. NGC 288 was not observed long enough to reach this luminosity limit of Lx ~ 4.0 x 10'" ergss"' in the 0.5-6.0 keV range. Its limiting 
luminosity is ~ 5.7 x 10'" ergs s"'. However, we keep this cluster's data, as a lower limit, since we have so few constraints on low-density 
clusters. 

c. No information about the optical counterparts to X-ray sources. 



a ve ry sparse cluster, which is evaporating its lowest-mass stars 
dHur lev et al. 2005). Open clusters don't survive very long. As 
a cluster evaporates its lowest-mass stars, it tends to retain its 
heaviest systems-binaries-which are more likely to be X-ray 
sources, b) M67 is a rather young cluster. Younger stars may 
produce more X-rays (• RandichligQTh . since they tend to be rotat- 
ing faster than older stars, c) A large fr action of binary syst ems 
are destroyed in globular clusters (see llvanova et al. I l2005h . m 
particular those with longer orbits. RS CVn systems involve gi- 
ants that are spun up by stellar companions. These systems must 
be relatively wide binaries, in order to avoid the giant swallow- 
ing its companion as it evolves; but such wide binaries are de- 
stroyed in globular clusters. Thus there are fewer RS CVn bina- 
ries in globulars. Since RS CVn binaries tend to be brighter than 
BY Dra binaries (main-sequence ABs), low-density clusters can 
have startlingly high X-ray luminosities per unit mass. 

To summarize the results of this paper, the number of X-ray 
faint sources with Lx,o.5-6 keV > 4 x 10^" ergs s"' found in M71 
is higher than the predicted value on the basis of either the col- 
lision frequency or the half mass. We suggest that those CVs 
and ABs in M71 are primordial in origin. The last interpreta- 
tion above may explain the X-ray overabundance of low-density 
clusters like M71, where fewer primordial binaries may have 
been destroyed through binary interactions. Study of other low- 
density globular clusters will help us to better understand their 
evolution and dynamics. 
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